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AbBtract 

This paper presents an analytical investigation 
of higher order mode propagation in a nonuniform 
circular duct without mean flow. An approximate 
wave equation is derived on the assumpLlons that the 
duct cross section varies slowly and that mode con- 
version is negligible. Exact closed form solutions 
are obtained for a particular class of converging- 
dlverging circular duct which is here referred to 
as " circular cosh duct ." Numerical results are 
presented in terras of the transmission loss for the 
various duct shapes and frequencies. The results 
are applicable to studies of multimodal propagation 
as well as single mode propagation. The results 
are also applicable to studies of sound radiation 
from certain types of contoured Inlet ducts, or of 
sound propagation in a converging-diverging duct of 
somewhat different shape from a cosh duct. 


List of Symbols 

a effective length of converging-diverging 

section of duct 

b^ duct radius at throat (x = 0) 

b(x) duct radius 

b_ constant radius of lefthand side uniform 

duct e’eraent 

b constant radius of rlghthand side uniform 

duct element 

c sound speed 

F hyporgeometrlc function 

J Bessel function of first kind of order m 

m 

k free space wave constant, o)/c 

k_ propagation constant of a mode in lefthand 

side uniform duct element, 

k propagation constant of a mode in righthand 

side uniform duct element, 

m circumferential mode number, or separation 

constant 


t time variable 

TL transmission loss (dB) 

a dimensionless parameter, see Eq. (19) 

X axial coordinate, see Fig, 1 

a separation constant or eigenvalue 

a eigenvalue, n-th zero of j'(a) 

Miu m 

0 contraction ratio, h_/h^ 

r Gamma fu.iction 

Y otitoff ratio of mode referenced to inlet 

kb_/a 

t dimensionless frequency parameter, 

kbpCl - p/Y) 

t| normalized variable, see Eqs. (9) and (9a) 

0 polar angle in local spherical coordinate 

system, see Fig. 3 

half-cone angle corresponding to a duct 
segment, see Fig. 3 

(1 asymmetry parameter, 1/4 In [ (p^ - 1/ t ^)/ 

(r2 - 1)] 

V converging-diverging section length paramo' 

ter, a/b_ 

1 dimensionless axial coordinate variable, 

see Eq. (19) 

p radial variable in cylindrical coordinate 

system 

a constant parameter, see Eq. (20) and the 

following paragraph 

T ratio between inlet and exit radii, b /b 

*T* 

0 function, see Eq. (17) 

tp azimuthal angle around duct axis 

T acoustic velocity potential 

CO angular frequency 

1. Introduction 


N normalization constant, see Eq. (22) 

n integer used as radial mode number 

q a dimensionless parameter, see Eq. (19) 

R power reflection coefficient 

r radial variable in spherical coordinate 

system 

T 

power transmission coefficient 


One of the unique features of the noise field 
produced by axial flow fans or compressors is Che 
dominfj.'Ce of spinning modes in the fan duct. An 
additional feature in aircraft applications is that 
the c'-' jfiB section of a fan duct varies along the 
duct axis. Consequently, in order to understand 
the overall acoustic characteristics of a fan duct 
system, it is imperative to examine the propagation 
of spinning modes (or higher order modes) in non- 
uniform ducts. 


"Based on Consulting Reports submitted to Fratt & 
Whitney Aircraft (November 25, 1973 and 
February 28, 1974). 

**Aerospace engineer; member AIAA. 

'''Frofessor of Physics and of Aeronautics and 
Astronautics. 


The higher order mode propagation in non- 
uniform ducts has been previously investigated by a 
number of authors. Tlie widely used methods are 
numerical methods, semi-numerical methods,^ vari- 
ous perturbation methods such as the WKB method or 
multiple scale variable methods^ and the Born ap- 


proxlmaClon.® One of the dtawlseka of numerical 
methods la the failure to provide compact expres" 
sions for the various quantities of interest and a 
corresponding lack of physical insight into the 
problem. Furthermore, in numerical methods the 
boundary conditions are satisfied only approximate- 
ly and the resulting error is not always easy to 
assess, perturbation methods have been used to 
find approximate solutions, but they still require 
extensive computation to yield numerical results 
for quantities of practical interest, such as the 
transmission coefficients. Furthermore, most ap- 
proximations used so far have been limited to first 
order perturbations with related limitations in the 
range of validity of the results, such a"' the limi- 
tation that the relative constriction be -t.mull. 

Also, it should bo mentioned that the fin»t order 
perturbations predict no attenuation above the cut- 
off frequency at the throat. Although on extension 
to higher order perturbations con be made formally, 
it is quite cumbersome to carry out, and it has 
rarely been done. 

In the present paper, we develop a new ap- 
proach with less limitation. The basic idea in- 
volved is to recognize that higher order mode propa- 
gation in a nonuniform duct is analogus to prob- 
lems in quantum mechanios dealing with matter wave 
propagation in a system with spatially varying 
potential energy. In the latter case, exact solu- 
tions ore known for certain barrier types of poten- 
tial.^ These solutions are used here to solve the 
acoustic problem for a class of converging-diverging 
circular ducts, which ore equivalent to the poten- 
tial barriers. Tlie duct shapes can bo varied by 
means of three independent duct parameters, and 
cover a wide range of ducts of ^'ractical interest. 
The final results include wave functions and trans- 
mission and reflection coefficients, all in a 
closed form. 


2, Circular Cosh Duet 


A typical nonuniform circular duct under con- 
sideration, which will be referred to os " circular 
cosh duct ." is composed of two asymptotically uni- 
form circular duct elements which are smoothly 
coupled through a converging-diverging section, as 
illustrated in Fig. 1. The shape is completely 
determined in terms of the radius, b(x), which is 
given by 

/ b . 2 

(b^) “ • =°sh(2n) + - De^l^ 

X j^coshV • sech^ “ sinh(2n) • tanh 

(1) 

where the parameters are defined in the symbol list. 
For a symmetric circular cosh duct (t = 1, |i = 0), 
this equation reduces to 

(bfe) = 1 (f). 

The duct shape given by Eq. (1) can be adjusted by 
means of the three dimensionless parameters p, v, 

T, allowing one to vary (i) the contraction ratio, 
p, (ii) the effective length of the converging- 


diverging section, V, and (Hi) the ratio between 
the inlet and exit radii, T. Equation (1) should 
cover a wide range of converging-diverging circular 
ducts of practical interest. Various circular cosh 
duct shapes with fixed p and v are displayed in 
Fig. 2, 


It should be mentioned here that the results 
of the present investigation can be used to study 
the acoustic characteristics of a variety of 
converging-diverging ducts with slowly varying 
cross section, which may slightly differ in shape 
from the circular cosh ducts. The supporting argu- 
ment for this statement is os follows: The fre- 

quency range where our understanding of acoustic 
properties of converging-diverging ducts is lock- 
ing, is near the cutoff frequency of a relevant 
duct mode. In this frequency range, the axial 
variation of wove phase is slow, and the propaga- 
tion characteristics are hardly affected by compara- 
tively small scale changes of the duct shape, as 
long as the cross section varies slowly. 


3. Wave Equation in a Non-uniform Circular Duet 

In this section, a horn equation will be de- 
rived for higher order mode propagation in a non- 
uniform circular duct under the following assump- 
tions: (t) Tlie slope, b', of the duct wall remains 

small, and (ii) The mode corivers ion (from one mode 
to another) is negligible. In addition, modes for 
such a non-uniform circular duct will be defined. 


Tire equation is derived first for a small seg- 
ment of the duct using locally suitable spherical 
coordinates . and is then transformed into a form 
which Involves the duct parameters and the axial 
coordinate x. Let us consider a small segment of 
the duct and local spherical coordinates (r, 0, q>) 
as illustrated in Fig, 3. The duct segment is so 
short that is may be regarded conical, that is, b' 
hardly changes within the duct segment. The origin 
of the spherical coordinate system is chosen some- 
where on the axis such that a coordinate surface 
0 •= 0Q (cone) tangentially contacts the wall of the 
duct segment. The half cone angle 8q, is equal to 
b' to second order approximation. 


With the assumption of harmonic time depen- 
dence e"^b the wove equation for spherical coor 
dlnates is 


_1_ 

,2 



+ 


-7^ ^ (sin 0 
sin 9 hQ \ 



, 1 

'00 
Bind 9cp 


+ k^T 


0 , 


( 2 ) 


where T is the acoustic velocity potential and 
k = uj/c. With the substitution 


T = H(r)Y(0)Q(qi), 

Equation (2) is separated as follows: 

v2-l 


A 2 

, 2 ^ r dr ^ 
dr 


k^H- 


fe) 


H = 0, 


sin 0 




(sin 0 + 
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( 4 ) 

(5) 
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AS + m^Q . 0. 
d<p 


CotcQppondlng to each combination o£ m and 
(S) “inn> '^1'“ Bolutlon 1* uniquely defined In terms of 
the product 


Here a and m arc the separation constantsi and (a q). 

0Q has been Inserted for convenience. 


We consider Eq. (6) first, the solution to 
uhleh Is 




e^"^. 


(7) 


Here m must be an Integer for Q to be a single 
valued function, and Is known as the circumferential 
mode number. It la obvious that this solution Is 
valid throughout the duct as long as the cylindri- 
cal symmetry Is maintained. 


Equation (5) Is an equation for the associated 
Legendre functions. These are not, however, con- 
venient for the present problem. The main diffi- 
culty arise? irv determining eigenvalues and in 
transforming thum from the local coordinates to the 
duct coordinates. Therefore, a simplified form of 
Eq. (5) Will be used here. In the region of inter- 
est 0 < S < 00, one Can replace sin 0 by 0 to 
the second order approximation in b'. Equation (5) 
can then be written as 


This eigenfunction Is defined throughout the duct, 
q being given by Eq. (9) for a conical segment or 
by Eq, (9a) for a cylindrical segment. Consequent- 
ly, a mode, which Is in one-to-one correspondence 
with an eigenfunction, has been defined for a non- 
uniform duct. 

We now return to Eq. (4). This equation, when 
transformed to the duct coordinates, Is the one 
which governs propagation of a mode In the duct. 

Tlie necessary transformation Is accomplished by re- 
placing Che local coordinate r with the duct 
parameters and the axial coordinate x. To this 
end, we use the following relations: 



(13) 

? ^ dj’ 

(14) 

~ A 

“ b" 

(15) 



( 8 ) 


These equations are valid for a point within the 
duct segment, and are accurate to the order of 

(b')L 


Tills Is the Bessel equation. 

It is convenient Co introduce a normalized 
variable q defined as 

n - -f, (9) 

0 

where 0 ^ q ^ 1. Note that, in the limit 0 q -• 0, 
the duct segment becomes cylindrical, and Eq. (9) 

Is replaced by 

n ' b’ 


where p la t%e radial coordinate In the cylindri- 
cal coordinate system. 

With the substitution of Eq. (9), Eq. (8) is 
transformed into 


d^Y 1 dY 

dq2 n dq 





Y = 


0 . 


( 10 ) 


On inserting Eqs. (13) to (15) into Eq. (4), 
we obtain 

This is the equation for higher order mode propaga- 
tion In a nonuniform circular duct. Its solutions 
will be sought In the following section. 


4. lIlRher Order Mode Propagation In a 
Circular Cosh Duct 


In this section, we seek solutions to Eq. (16) 
for the circular cosh ducts described in Section 2. 

It is convenient to Introduce a new function 
$(x), which is defined by 

H(x) = (17 


Note that q Is not a local coordinate variable, 
and unlike Eq. (8), Eq. (10) is valid throughout 
the duct as long as the assume Mons (1) and (11) re- 
main valid. The physic, illy acceptable solution to 
Eq. (10) is the Bessel function of the first kind. 
Thus we have 

Y = J^(o,q). (11) 

From the boundary condition at the duct wall, eigen- 
values are determined as follows 

a ” ~ 8,1,2 . . . (12) 

being Che n-th zero of dm (“■) • 


On inserting this into Eq. (16), we obtain 



where bb'7a has been discarded on the assumption 
that it is negligible compared with unity. It can 
be readily shown that this equation is also valid 
for the fundamental mode propagation (a = 0) . 

Equation (18) is similar, in form, to the one 
dimensional Schroedinger equation, with k^ and 
(a/b)2 corresponding to the total energy and the 
potential energy respectively. 7 With the substitu- 
tion of Eq. (1), Eq. (18) is written as 


3 


di 


0 r 2 2 2 

■j + I q - Vjj(coah (i • sech t 


- sinh(2)i) • tonli J)j 0 » 0, 


(19) 


where 


X - ua 


q^ >• (ka)^ + Vjj cosh(2(i) - (apv)^, 

V„ - T 1) . (c..v)^ 

Equation (19) can be aolved exactly in accor- 
dance with the procedure given in Section 12.3 of 
Uof. 7. With the incident wave coming from 
X " the Bolution is given by 

I (k, - k_)o5 I (k + kja 

0 - N e‘‘ ^ (2 cosh 6)^ ^ 


Yi . i 

\2 2 


1 i 


K F I j j (k_^a + k_a - J ~ 2 (k^*" + k_a + 0), 


1 - Ik 


e^^'q. i)‘ 


(20) 


Here F is the hypergeometric function, N a con- 
stant to bo dotortnlned, and 


a “ cosh^n - 1. 


This solution, 6 being replaced by x, can be 
written asymptotically at x « ±co as follows 

4 ^ N e’^k+)xa „ik^.x 


(21a) 
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The solution in Eq. (20), Eq. (17) and the 
corresponding eigenfunction will give ,ne complete 
solution in the circular cosh duct; For the inci- 
dent wave given by 




isq) ^ j e^A.’i at X •• -w, 
m' mn " ’ ’ 


(23) 


the resultant wove in the duct is 
im9 


m mn ' b 


+ e 


1 f (k+ - kj (X - pa) 
b ® 

- (f - *)■ 


"* f (k^ + k_)a 


X F || - I (k^o + k_a - 0), I - I (k_^a + k_a + a), 


1 - ik_j^a; 


(f--) 


(24) 


+ ly 


The reflection and transmission coefficients 
for the wave a>'.plltude can be readily computed from 
Eqs. (21a and b^ In the case of no mean flow, the 
reflection and transmission coefficients for the 
acoustic power are obtained as the absolute squares 
of the respective amplitude coefficients. In the 
following diseussiens, we will use the acoustic 
power reflection and transmission coefficients, 
which are given by 


T = 


cosh['n'a(k, - k )] + cosh(TrO) 

*r *'** 

coshtTra(k, + k )J + cosh(TrO) ’ 

T 

2 slnh(Trk^a) • sinh(Trk_a) 
cosh[Tra(k^ + k_)J + cosh(Tro')' 


(25) 


( 26 ) 


In deriving these equations, we have used the 
following properties of the P function. 


, N • P(1 - Ik, a) • r(-ik a) • e^A.x 

' *r *** 

r (I - I (k_^a + k_a + 0 )) • P (I - i (k^a + k_a - o)) 
N • P(1 - ika) • P(ik a) ■ e^kjxa g-lk_x 

*r *“* 

r (I - I (k^a - k_a + 0 )) • P (i - i (k^a - k_a - o))’ 




ttz 


'(iH-iz)f 


cosh (ttz) ' 


Note that the results given in Eqs. (25) and (26) 
satisfy energy conservation, that is 


(21b) 


R + T = 1 


(27) 


where P is the gamma function. The asymptotic 
solution in Eq. (21a) corresponds to the trans- 
mitted wave, and the first and the second terms of 
Eq, (21b) correspond to the incident and the re- 
flected waves respectively. Requiring the incident 
wave to have unit amplitude, we have 

b_r (I - I (k^a + k_a + g)) T (I - I (k^a ■)- k_a - a)) 

P(l-ik a) -P(-ik a) 

*r “ 

( 22 ) 


5. Numerical Results and Discussions 

Some numerical results are presented in 
Figs. 4 to 7. Displayed in the figures is the 
acoustic power transmission loss (TL) for various 
values of the duct parameters and a frequency 
parameter (y or 5)- I^ke TL is defined as 

TL(dB) = -10 logj^Q(T). (28) 

In Fig. 4, the TL of the (1,0) mode is 
plotted as a function of y, the cutoff ratio of 
the mode at the. duct inlet, defined as 
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kb 

Y - — (29) 

Note that the value of y la not leaa than unity 
for propagating Incident wave mode. The duct 
parameters Included In the figure are various val- 
ues of 0 and V » T *• 1. As expected, the Th de- 
creases with Increasing frequency, In a similar way 
for different values of 0. Ttte rate of the TL 
change (dTl/dy) Is the greatest In the vicinity of 
V " 01 which corresponds to the cutoff frequency of 
the mode at the throat. Note that, at y ■* 0, the 
TL Is nonzero with a value of about 2.5 to 3 dB. 

This result Is In contrast to that of the first 
order perturbation solutions, which predict no 
attenuation for y > 0.9>0 

In Klgs. 5 and 6, the TL of various modes la 
plotted as a function of for various values of 
the duct parameters. Here f, is a frequency 
parameter which Is defined as 

5 - kb^ (l - f)- (30) 

Note that 0/y Is the Inverse of the cutoff ratio 
at the throat, and that f, “ 0 corresponds to the 
cutoff frequency at the throat. 

Each curve (narrow stripe) in these figures 
Includes the TL of ten randomly selected modes rang- 
ing from the (2,0) mode to the (8,5) mode. This 
shows that the f,- dependence of the TL is almost 
the same for all the higher order modes. This ro= 
suit is remarkable, and is useful especially for 
studies of multimodal propagation In a converging- 
diverging duct. In some cases, one may not need to 
Identify individual modes, An example of such 
cases can be found In conjunction with a study of 
multimodal radiation from a uniform duct.® The 
acoustic power distribution vs. y in the duct can 
often be inferred from the measured radiation pat- 
tern. For a known ydlstribution of the acoustic 
power, the present result can be immediately ap- 
plied to determine the effect of the converging- 
diverging section on the multimodal propagation. 
Recall that the parameter ^ contains y explicit- 
ly, not the eigenvalues of modes. 

Figure 5 includes calculations for 0 = 1.05, 
1.25, 1.5 with V “ T * 1. Tb-: result shows that, 
for 5 < 0, the TL is more sensitive to the fre- 
quency for the smaller ’'alue of 0. We also notice 
that the calculations for different modes collapse 
better as 0 increases. 

Figure 6 Includes calculations for different 
values of v and t, with 0 = 1.25. We first com- 
pare the curves (narrow stripes) for (v = 1, T = 1) 
and for (v = 1, T = 1.5) in order to see the effect 
of T. For ^ < -0.15, the TL for T = 1.5 is 
smaller than that for T *= 1. This result is due 
to the fact that the duct segment In which the modes 
are cutoff, is shorter for T = 1.5 than that for 
T = 1. (See Fig. 2). However, the two curves do 
not show much difference for ^ > -0.15. This 
point will be further discussed later. 

We now compare the curves for (v = 1, T “ 1.5) 
and for (v = 1.5, 1.5) in Fig. 6 in order to 

see the effect of v. For ^ < 0, the TL is larger 

for V = 1.5 than that for v = 1. This result is 

the direct consequence of the fact that the larger 
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values of V represents the longer converging- 
diverging section. The length of the duct segment 
In which the modes are cutoff, Increases with In- 
creasing value of V. 

In Fig. 7, the TL of the (1,0) mode la plotted 
as a function of T for the varioua values of y, 
V, and 0, For y > 0, the TL hardly depends on t. 
For y < 0, Che TL slightly decreases first with 
Increasing values of t, and asymptotically reaches 
a constant value near t “ 1.5. It remains un- 
changed for further Increase of the value of x* 
This finding suggests that the results of the pre- 
sent Investigation can be used for studies of high- 
er order mode transmission from a contoured Inlet 
duct. 


TO account for the last statement, we consider 
a contoured Inlet duct as illustrated In Fig. 9. 

The Inlet duct is produced from a circular cosh 
duct with a large value of x. The circular cosh 
duct is terminated at a distance where the duct 
radius is 1.4b_. Tlic cutoff frequency at the 
termination, is smaller than that of the lefthand 
side uniform duct element by the factor of 1.4. 

Tlius, for a propagating incident wave mode (y > 1), 
the reflection due to the duct termination is neg- 
ligible.® It follows then that the TL for the con- 
toured inlet duct is expected to be the same as 
that for the full circular cosh duct. 

It should probably be mentioned that the slope 
b ' can be fairly large for a large value of x In 
a portion of the diverging section. Although it 
seems to violate the assumptions used in the analy- 
sis, the numoricol results remain valid. The sup- 
porting argument is that, in the first place, the 
slope remains small within some distance downstream 
from the turning point where the mode changes from 
cutoff to cuton. In the second place, the sound 
attenuation takes place mostly in the duct segment 
‘here the mode is cutoff, whereas the diverging 
section, like a loud-speaker horn,!® helps the un- 
attenuatod portion of the sound to bo radiated. 

The radiation efficiency is not sensitive to the 
slope ns long as the duct divergence is smooth and 
the frequency is not close to the cutoff frequency 
at the termination.® Lastly, oven a fairly largo 
slope, e.g., b' <= 1, is not so large as to invali- 
date the approximations made in the analysis. For 
instance, for b' = 1, = Tr/4 and sin Sq ° 0.707; 

thus, the error involved in replacing sin Sq 
9o is about 11%. 


6 . Concluding Remarks 

In an attempt to Improve the understanding of 
the acoustic characteristics of a fan duct system, 
we have investigated higher order mode propagation 
in a particular class of converging-diverging cir- 
cular duct, called "the circular cosh duct ." The 
duct shape can be adjusted by means of throe duct 
parameters, covering a wide range of converging- 
diverging ducts of practical interest. 

An approximate wave equation has been derived 
and exact closed form solutions have been obtained. 
No mean flow effects have been included. The ex- 
pressions for the reflection and transmission coef- 
ficients of a mode are simple. With an appropriate 
choice (^) of frequency parameter, the results have 
been shown to be nearly independent of individual 
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modes, and can, therefore, be immediately used for 
studies of multimodal propagation. The results are 
also applicable to studies of sound radiation from 
certain types of contoured inlet ducts, or of sound 
propagation in a converging-diverging duct which 
differs somewhat from a cosh duct in shape. 
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Figure 4. - Pwer transmission loss of (1,0) mode versus frequency 
parameter for various values of p, with u » t - 1, 


Figures. - Transmission loss versus for p*1.05, 
1.25, and 1.5, for u » t ■ 1. Each curve (stripe) in- 
cludes 10 modes ranging from (2.0) mode to (8,5) mode. 


(U'1,T •=1.5)-' 


(v = 1. 5, T = 1. 5) 


■{V = 1, T ■= 1) 


Figured. - Transmission loss versus for v = l, t' 
V = 1, T = 1. 5; V ■= 1. 5, T = 1. 5 with p =• 1. 25. Each 
curve (stripe) includes 10 modes ranging (2,0) mode to 
(8, 5) mode. 
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16. Abstract 

This paper presents an analytical Investigation of higher order mode propagation In a nonuniform 
circular duct without mean flow. An up}, roxlmate wave equation Is derived on the assumptions 
that the duct cross section varies slov/iy and tliat mode conversion is negligible. Exact closed 
form solutions are obtained for a particular class of converging-diverging circular duct which 
Is here referred to as " circular cosh duct. " Numerical results are presented In terms of the 
transmission loss for the various duct shapes and frequencies. The results are applicable to 
studies of multimodal propagation as well as single mode propagation, The results are also 
applicable to studies of sound radiation from certain types of contoured Inlet ducts, or of sound 
propagation In a converging-diverging duct of somewhat different shape from a cosh duct, 
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